Brownian motion of the electron and the Lamb shift at finite temperature 
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By enhancing electron position fluctuations, equilibrium electromagnetic radiation modifies the 
potential for an electron in a Hydrogen atom. This can have significant effects for weakly bound 
states and especially at finite temperature. This implies a 2% correction to Bethe's value for 2Si/2 — 
2P1/2 Lamb shift for weak fluctuations, but the effect is an order of magnitude larger for strong 
fluctuations where it provides direct measure of the proton diameter. 



PACS numbers: 31.30.jf, 05.40.-a, 05.40.Jc, 03.65.Sq 

In 1948, following Bethe's classic calculation [l] of the 
experimentally observed energy splitting of the 2Si/2 and 
2Pi/2 levels of Hydrogen, the celebrated Lamb shift Q], 
Welton J3:| presented an elegant semiclassical explanation 
of the effect in terms of fluctuations of the electron posi- 
tion superimposed on its orbital motion. The Brownian 
motion of the electron is due to the coupling to zero- 
point fluctuations of the electromagnetic field. Welton 
observed that this generates an effective repulsive barrier 
at the origin which he treated via first-order perturbation 
theory. The goal of this paper is to point out that the 
assumption that this induced potential is weak is not sat- 
isfied as well as one might expect. This leads to a small 
correction to the Bethe value of the 25i/2 ~ 2Pi/2 level 
splitting. The effect is an order of magnitude larger for 
strong electron fiuctuations which will be the case for a 
sufficiently weakly bound electron, and especially at finite 
temperature. Since these findings have non-relativistic 
origin and Welton's analysis, while reproducing the bulk 
of the Lamb shift, is unaffected by the many-particle and 
relativistic effects jj], the full quantum electrodynamics 
treatment [5| is unnecessary. 

Fluctuations of the electron position induce changes 
both in the kinetic and potential energy. The former can 
be ignored because they are the same for both the S and 
P states; only the change in the potential energy U{r) 
needs to be considered. When the electron trajectory 
is perturbed so that r — ^ r + (5r where Sr is an isotropic 
fiuctuation with zero average < 6r >— and finite mean- 
square displacement < 6r^ >, the potential energy felt 
by the electron will be given by [3| 

<U{r + Sr)>^exp(^ <5r^ >\/Au{r) (1) 

If the proton would be a point charge e, then C/(r) — 
— e^/r and the extra induced potential felt by the electron 
would be [3| 



SU{r) ~ — e^ < Sr^ > S{r) 



(2) 



Such three-dimensional delta-function potential, in con- 
trast to its one-dimensional counterpart, has a formal 
meaning and can only be viewed as a perturbation and 



only to first order. The true potential at small distances 
is a barrier of width tq, the proton radius, and height 
6U ~ e^ < (5r^ > /rg; whether this barrier can be re- 
garded as weak or strong requires separate investigation. 
In order to do so properly, we will model the proton by a 
uniformly-charged ball of radius tq. Then for r > ro the 
potential energy of the electron is U{r) = — e^/r while 
for r ^ rp it is given by 

C/(r) 



2ro 



^0 



(3) 



As a result the electron experiences an attraction that 
is weaker than —e^jr at small distances, which increases 
its energy. The corresponding energy shift is quadratic 
in the proton radius rg [6j regardless of the model of 
the charge distribution inside the proton. Since r^/aB = 
1.66 X 10~^ [3], where qb = T? jrae^ is the Bohr radius, 
the finite size correction to the energy due to ([3]) will be 
neglected (see below) compared to the fluctuation effect. 
Substituting ([3]) into ([T]) we observe that due to the 
position fluctuations of the electron the potential energy 
acquires an extra contribution given exactly by 



m{v) 






(4) 



which is a constant height barrier operating at r ^ ro. 
The strength of such a barrier can be judged by forming 
the dimensionless combination 

1/2 

x = \ r-z:''' = (— i;^-— 1 (5) 



ft2 



ro 



The consequence is that because as/rg :3> 1, the condi- 
tion of weak electron fluctuations < Sr'^ > /a% ^ 1 alone 
does not guarantee the weakness of the induced barrier 
Q. The crossover between the regimes of weak A ^ 1 
and strong A ^ 1 barriers occurs at < Sr >^ /a\ ~ 
Ta/cLB ~ 10~^, i.e. when position fluctuations are small. 
We thus face a Coulomb problem modifled at small dis- 
tances by the presence of the barrier @ ; principally the 
S states are affected by the barrier. The upper part of the 
5-state energy spectrum of a generic modifled Coulomb 
problem is described by the Rydberg formula 

^„(/i) = const - , _ ,^ {a.u.) (6) 



where we employed atomic units (a.u.), the const is com- 
mon to the S and P states and n is the quantum defect 
[8|. The latter represents convenient dimensionless mea- 
sure of experimental splitting between the S and P levels 
determined by the difference of © evaluated at finite and 
zero /i. The quantum defect can be expressed in terms 
of the S'-wave scattering length a^ of the r ^ vq poten- 
tial only [8|. In the present case of the constant height 
potential barrier (JH) the scattering length is given by 



ro 1- 



tanh A 



A 



(7) 



Since a^ ^ tq for all A and tq -C a^ , the scattering length 
also satisfies the condition Og <^ as (which is also true 
for any reasonable model of charge distribution inside 
the proton). This is the range of applicability of the 
perturbation theory in as /as, where the quantum defect 
is given by [3, [3] 



A* = 



2qs 
as 



(8) 



In the limit of a weak barrier A ^ 1 (with A given by 
Eq.©) the scattering length ^ approaches a limit that 



is independent of the proton radius tq : 



roA73 =< 



(5r > /ZaB <IC rg. The corresponding quantum defect 



2<5r^ > 



Mo 



(9) 



directly measures electron's mean-square fluctuation. 
This conclusion is equivalent to the result of Welton [3i] 
that follows from the |/i|/n ^ 1 expansion of the Ryd- 
berg formula (O ; this is the regime of applicability of the 
perturbation theory in the induced potential ([2])- 

In the large-barrier limit A 3> 1 the scattering length Os 
([7]) approaches the size of the proton r^ with the quantum 
defect 



A*c 



2ro 
as 



-3.32 X 10"^ 



(10) 



directly measuring the negative of the proton diameter. 
We note that the weak ([9]) and strong pO)) barrier results 
are insensitive to the model of the charge distribution 
inside the proton. 

We now proceed to a calculation of the mean-square 
fluctuation < 5r^ > which will tell us what regime the 
real system may belong to. Consider a non-relativistic 
(|r| <^ c) harmonically bound electron of mass m driven 
by a fluctuating time-dependent electric field E(i). The 
force on the electron corresponding to this field is f = 
— eE and the equation of motion for the electron - in- 
cluding radiation damping - is [lOl] 



mv 



2e2. 
3c3' 



mujf^Y 



f(t) = -eE(t) 



(11) 



where wq is an oscillation frequency representing the fre- 
quency of motion in the Bohr orbit. This equation differs 
from that employed by Welton in one crucial respect: 

As written, Welton's equation of motion of a free 
(wq = 0) electron mv = — eE(i) with E(i) represent- 
ing the fluctuating equilibrium electric field contradicts 
the fluctuation-dissipation theorem (FDT) [ll|, |l2| since 
it contains fluctuations but no dissipation. It was indeed 
later recognized by Callen and Welton (llf that applica- 
tion of the FDT to Eq. ([TT|) leads to the Planck formula 
for the spectral energy density of the radiation. Thus the 
radiation damping term — (2e^/3c3)r in the equation of 
motion (jlip is both appropriate and necessary. 

Introducing a spectral decomposition of the fluctuat- 
ing quantities according to the conventions of Ref. jl^] , 
r(t) = /^ Ti^e"*'^*!^, the equation of motion (ITT]) can 
be written in the spectral form as r^^ — A(uj)i^ where 
A{uj) = A'{uj) + iA"{uj) is the generalized susceptibility: 



1 



A{u:) 



m{ujQ — Lo^) — i 



2,, ,3 



2e^uj 
3c3 



(12) 



This equation is applicable for frequencies that are sig- 
nificantly smaller than the inverse time of travel of light 
through a distance of the order of the classical radius of 
the electron, ui <C nu? j e^\ ignoring this constraint of in- 
ternal consistency of classical electrodynamics is known 
to lead to physically absurd results [10| . For our problem 
this condition is practically irrelevant because quantum 
effects come into play at frequencies of the order rac? jh 
137 times smaller than the clas- 



hcjt 



which is a~^ 
sical limit [lO|. 

The generalized susceptibility p^ determines, via the 
FDT, the mean-square displacement of the electron jl2| : 



< 8v^ >= 



— / ^"(c^)coth— dw 

TT Jo 2J 



(13) 



where the factor of 3 accounts for the three components 
of the displacement vector r and T is a temperature. In 
the neutral limit e — > the integral (|13p can be evaluated 
with the result 



< 8x' >= 



3h hojo 

coth 

2mujQ 2T 



OT-l \ J 



that can be recognized as Bloch's formula for the mean- 
square displacement of an undamped oscillator 12 1. 

For a free electron ujq = 0, and at T = Eqs.([T2]) and 
(HSl) lead to the result 



< Sv' >= 



2a-^ 



did 



t^ + (2eV3mc3)2a;3 



(15) 



This expression differs from Welton's only in the presence 
of the (2e^/3mc3)2a;3 term in the denominator. The fre- 
quency integral is logarithmically divergent at the lower 
limit; for a bound electron the effective infrared cutoff 



is Wo 3]. In contrast to Welton's result, the integral 
in (|15p is ultraviolet convergent; the effective frequency 
cutoff is set by mc'^/e^, the range of applicability of clas- 
sical electrodynamics. Welton recognized the necessity 
for a physical ultraviolet cutoff which he set at a fre- 
quency of the order mt? jh where quantum-mechanical 
effects come into play. Therefore it is appropriate to ig- 
nore the (2e^ lZm(?)'^ufi term in the denominator in (fT5|) 
and instead set the upper integration cutoff at rru? jh. 
This leads to Welton's result for the mean-square dis- 
placement of a bound electron 



< Jr^ > 2o? 



duj 



2o? mc 

■K hujQ 



2a^ , n^ 
In^; 



TT a 

(16) 
which has logarithmic accuracy. Comparing the ojq- 
dependence here and in Eq. (fT4| at T = we conclude 
that the fluctuations of charged oscillator are dramati- 
cally suppressed. In the last step in (J16p we employed 



UJQ 



hasn^ 



(17) 



(Bohr's correspondence principle) for the frequency of 
motion in Bohr's n-th orbit applicable for n 3> 1. For 
the 2S state and reasonable cutoffs the logarithm in p^ 
is approximately 8 [l| and then Eq.(fT6| tells us that at 
zero temperature the root-mean-square fluctuation of the 
electron in the Bohr orbit is about 700 times smaller than 
the Bohr radius. Even though the mean-square displace- 
ment grows with the principal quantum number, the n- 
dependence is weak: for n — 50 the root-mean-square 
fluctuation is only about 1.5 times larger than for n = 2. 
For the same 25' state substitution of p^ into © leads 
to the quantum defect /io = — 1.32 x 10~^ whose mag- 
nitude is 25 times smaller than the strong barrier result 
(fT0| . This however assumed that the induced barrier can 
be regarded as weak. In order to verify whether this is 
the case we compute the effective strength parameter ([5]) 
with the result A w 0.35 which is smaller than unity (but 
not by a lot). Then a recalculation of the quantum de- 
fect ([8]) with the scattering length evaluated according 
to Eq.© with A = 0.35 gives fi = -1.29 x 10"^ which 
is 2.27% different from Bethe's value. This 2.27% differ- 
ence has the same order of magnitude as the contribu- 
tion into the Lamb shift due to vacuum polarization [5[ 
and is two orders of magnitude larger than the leading 
flnite-size correction contributing an amount of the or- 
der — (ro/as)^ [y] into the quantum defect. We also note 
that what we believe is the correct non-relativistic 2S 
quantum defect {fi = —1.29 x 10~^) depends only weakly 
on the assumption of a uniform charge distribution inside 
the proton. 

It would be important to verify whether this effect 
could resolve the ongoing controversy regarding the size 
of the proton because presently accepted value of the pro- 
ton radius [7] is extracted from precision spectroscopy of 



atomic Hydrogen. However supposedly more accurate 
measurement involving muonic Hydrogen gives approxi- 
mately a 4% smaller value of the radius [13| . Our analysis 
is not applicable to the case of muonic Hydrogen where 
the bulk of the effect is due to polarization of electron 
vacuum [5]: the muon Bohr radius 0^/207 is within the 
as/137 range where polarization of the electron vacuum 
substantially modifies the Coulomb law. Thus the muon 
experiences an attraction that is stronger than — e^/r at 
small distances which decreases its energy. 

It is also important to ask if it is possible to observe 
much larger deviations from the Bethe- Welton value ([9]) , 
specifically the impenetrable barrier limit (|10l) . A gener- 
alization of our analysis to the case of a Hydrogenic ion 
of atomic number Z <i; 137 and atomic mass A predicts 
that with logarithmic accuracy this amounts to multipli- 
cation of Eq.dH) by Z^/^A'^/^. For sodium ion (Z = 11, 
A — 23) this increases the effective strength parameter 
([5]) by a factor of 2. However, significantly larger values 
of A can be realized at finite temperature and n ^ 1 due 
to enhancement of the electron fluctuations. 

For a free electron (wn = 0), a counterpart to Eq. p^ 
was given by Moore [IJ] who concluded that thermal 
fluctuations have negligible effect on the 25*1/2 ~ '^Pi/2 
Lamb shift. Although we agree with such an assessment, 
Moore's account of the electron binding was incorrect; 
we also disagree with his choice for the ultraviolet cutoff. 

At flnite temperature the mean-square displacement 
(fT3| can be approximately computed by setting coth x w 
1/x for X < 1 and coth a; w 1 for a; ^ 1, which sepa- 
rates the fluctuations into the classical (a; <C T/h) and 
quantum (a; ;» T/h) ranges |12J : 



< Sr' >. 



6T 



A"{lo) 



did 



m 



2T 
h 



A"{uj)duj 



(18) 

Assuming that the temperature is high enough [hijQ <C 
T) allows us to calculate both integrals: 

In the flrst integral the upper limit can be re- 
placed with inflnity and the integral, with the help of 
the Kramers-Kronig relationship [12[, is evaluated to 
{tt/2)A'{0) = (7r/2)^(0) = TT/{2mujl). For the second 
integral we can neglect wq relative to lo in ()12p . so that 
A"{ijj) = 2e^ / im? c^ uj . Then the mean-square displace- 
ment of the electron will be given by 



< 5r^ > 



3T 



3Tn^ 



2a^ 

TT 

2o? 

TT 



In 



mc 
~T~ 



In- 



Ta2 



(a.u.) (19) 



where in the second representation we employed Eq. (|17p . 
The second terms in (fT9|) give the contribution of the 
quantum fluctuations while the first terms can be rec- 
ognized as the classical mean-square displacement of the 
position of a three-dimensional harmonic oscillator (same 
as the T ^ hcuQ limit of the Bloch formula ([T?1) V 



Comparing Eas. p6|) and ((19)) we note a stronger di- 
vergence of the mean-square displacement in the latter 
case as the free-electron limit wq — s> (or n — > cx)) is 
approached. It is clear that the mean-square displace- 
ment cannot grow with n without bound as predicted 
by Ea.([T^. Indeed, the concept of the orbit retains its 
meaning only as long as the root-mean-square displace- 
ment is much smaller than the orbital size n^{a.u.). This 
leads to the constraint Tn^ ^ l(a.u). It seems plausible 
that the states satisfying Tri^ > l{a.u.) will be destroyed 
by thermal fluctuations; this means that at room temper- 
ature T « lO^^a.u. the n > 30 states are unobservable. 

Moreover, one has to make sure that during an ex- 
periment the state in question does not decay, which 
requires that its lifetime 15| n'^/{a^T){a.u.) is signifi- 



cantly longer than the Kepler period n^{a.u), thus im- 
plying n ^ a^^T~^{a.u.). However this condition is in 
practice irrelevant since it cannot compete for realistic 
temperatures with the requirement that the orbit is well- 
defined (Tn^ < l(a.u.)). Combining the Tn^ < l{a.u.) 
constraint with that of high-temperature limit T ^ huQ 
with ujQ given by (|17|) provides us with a range of princi- 
pal quantum numbers to which Eq. (|19p is applicable 



T 



-1/3 



<7i<cr 



-1/2 



(a.u.) 



(20) 



At room temperature this gives a range of n between 
10 and 30; at 40 K the range is between 30 and 100. 
With these values the mean-square displacement of the 
electron is dominated by the 3Tn^ term of P^ . 

The strong-barrier result pO)) should be observable in 
the regime of large electron fluctuations as described by 
Egs. p^ and ([^0]) . In practice we expect it to be valid 
in a much wider range of parameters because, as was 
already mentioned, the crossover between the regimes of 
weak and strong barriers occurs when the mean-square 
position fluctuation < 5r >^ /a^ « 10~^ is rather small. 

The results described above have their origin in the 
divergent behavior of the mean-square fluctuation of the 
free electron. Deeper understanding of what that means 
can be achieved by consideration of the correlation be- 
tween displacements of a free (wq — 0) electron at differ- 
ent moments of time, which are measured by the correla- 
tion function eft )=< (r(t) — r(0)) >, again determined 
by theFDT ^]^: 



C[t) = — / A"{uj) (1 - coswt)coth— dw (21) 
Ti" Jo 2T 

which at zero temperature and in the long time limit 
\n{mc^t/h) 3> 1 can be evaluated as 



C(t)«i^a|lnl^ 



(22) 



This implies that the electron executes a very slow sub- 
diffusive Brownian motion. 



At finite temperature the correlation function C{t) can 
be computed in the Tt/h ^> 1 limit with the result 



Cit) 



ATe' 



rt = 6Dt, 



D = 



2Te' 



(23) 



This can be recognized as the standard diffusion motion 
with diffusion constant D whose physical origin can be 
understood by observing that the result (j23p is purely 
classical. Therefore the Einstein relation between the 
diffusion constant D and static mobility B'{Q) must hold 
thus implying a very small mobility B'{0) = D/T — 
2e^/(3TO^c"^). Since the frequency-dependent mobility 
B{ijo) — B'{uj) + iB"{uj) is the coefficient of propor- 
tionality between the spectral components of the elec- 
tron's terminal velocity and the external force, we have 
B'{uj) = ujA"{uj). Employing Eq.dUl) with ujq = 0, and 
taking the static w — )■ limit we again find B'{0) = 
D/T — 2e^/(3m^c^). Therefore we conclude that the 
free electron mobility and associated diffusion constant 
have their origin in the radiation damping. 
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